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Abstract 

When two drops of radius R touch, surface tension drives an initially sin- 
gular motion which joins them into a bigger drop with smaller surface area. 
This motion is always viscously dominated at early times. We focus on the 
early-time behavior of the radius r m of the small bridge between the two 
drops. The flow is driven by a highly curved meniscus of length 2-7rr m and 
width A <C r m around the bridge, from which we conclude that the leading- 
order problem is asymptotically equivalent to its two-dimensional counter- 
part. An exact two-dimensional solution for the case of inviscid surround- 
ings [Hopper, J. Fluid Mech. 213, 349 (1990)] shows that A oc r m and 
r m ~ (^t/ 71 " 7 ?) hi [£7/(7?-??)] ; and thus the same is true in three dimensions. 
The case of coalescence with an external viscous fluid is also studied in detail 
both analytically and numerically. A significantly different structure is found 
in which the outer fluid forms a toroidal bubble of radius A cx rm at the 
meniscus and r m ~ (t"f / inn) In \t^/{n R)]. This basic difference is due to the 
presence of the outer fluid viscosity, however small. With lengths scaled by 
R a full description of the asymptotic flow for r m (t) <C 1 involves matching 
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of lengthscales of order ; ?*m , T~mi 1 and probably r m . 
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I. INTRODUCTION 



Considerable interest has been devoted recently to the breakup of free-surface flows into 
drops under the action of surface tension (Rallison 1984; Stone 1994; Eggers 1997). Here we 
investigate the complementary problem of surface-tension-driven coalescence of two drops, 
which is of fundamental importance in understanding the possible topological transitions 
in three-dimensional free-surface flows. For example, numerical implementations of merging 
(LaFaurie, Nardone, Scardovelli, Zaleski & Zanetti 1994) are based on phenomenological pre- 
scriptions for joining the two surfaces, without a fundamental understanding or description 
of the dynamics. 

Traditional applications of coalescence ideas include the description of two-phase disper- 
sions. As an important example, we mention phase separation in two-phase flows (Niko- 
layev, Beysens & Guenoun 1996; Bonnecaze, Martula & Lloyd 1998), where the velocity 
field induced by the merging of two drops entrains other drops, thus enhancing the rate of 
coalescence. Another classical problem connected with drop coalescence is sintering, i.e. the 
merging of a powder into a homogeneous material by heating. In many cases, in particular 
that of ceramics or glasses, bulk fluid motion is the dominant mechanism for coalescence 
and the dynamical process is known as viscous sintering. In a classical paper, Frenkel (1945) 
posed the problem of the merging of two spheres by slow fluid motion as the first step towards 
understanding the properties of the material that results from sintering. For a different case 
in which surface diffusion is the dominant mechanism of mass transport, the asymptotics 
of coalescence has recently been worked out (Eggers 1998), but such surface-dominated 
transport is very different from the bulk fluid motion of interest here. 

Much of the experimental and numerical work on coalescence in viscous systems is mo- 
tivated by the viscous sintering problem. An exception is an experimental paper (Bradley 
& Stow 1978) on the coalescence of water drops, but the low viscosity of water makes the 
motion very rapid and difficult to observe. On the other hand, by using a very high viscos- 
ity fluid, the motion can be slowed down as much as desired (Brinker & Scherer 1990) and 
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the experimental results agree very well with numerical simulations of the Stokes equations 
(Martinez-Herrera & Derby 1995). The only theoretical analysis of three-dimensional coa- 
lescence is the qualitative work by Frenkel (1945). Analytical solutions for two-dimensional 
coalescence (i.e. of parallel cylinders) have been obtained using complex variable techniques 
for the special case where the outer fluid is perfectly inviscid or absent (Hopper 1990, 1992, 
1993a,b; Richardson 1992). We show below that the three-dimensional problem has the 
same asymptotic behaviour as this two-dimensional solution at early times. Our main aim, 
however, is to address the more general case of coalescence with a viscous outer fluid, for 
which we find that the structure of the solution near coalescence is quite different from the 
case of an inviscid exterior, though there is again a parallel between the two-dimensional 
and three-dimensional problems. 

Part of the challenge in treating three-dimensional coalescence arises from the fact that 
it starts from a singular initial condition, shown in figure |l[ We assume that the drops 
are initially spherical, which is based on an underlying assumption of negligible velocity 
of approach and hence negligible hydrodynamic deformation before contact. We imagine 
that two such drops have just been joined along their symmetry axis by some microscopic 
mechanism to form a tiny bridge of radius r m . Evidently, the "meniscus" around the bridge 
will be a region of very high curvature, which drives the increase of r m with time. Our 
main concerns will be the time dependence of r m for very early times, and the shape of the 
interface and the flow field near the meniscus. We note that the ratio of the coefficient of 
surface tension 7 and the viscosity 77 gives a fixed velocity scale 7/77, and thus the expected 
Reynolds number Re = fryr m /r] 2 will be be arbitrarily small as r m — > 0, and the flow will 
initially be described by the Stokes equations regardless of the material parameters. 

In the next section, we set up an integral representation of the Stokes flow and split it 
into two parts: an outer region far from the meniscus, in which the shape is still close to the 
initial spherical condition, and an inner region near the meniscus in which the shape evolves 
rapidly. The dominant contribution to the velocity field comes from the high curvature of 
the meniscus, and its amplitude is determined by the lengthscale A of this curvature. The 
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main task is thus to find the structure and scale of this inner solution. Since the inner 
solution is determined by the local curvature and not the global shape, the results obtained 
here are not restricted to the simple spherical geometry shown in figure []]. In section III we 
show that only a negligible fraction of the fluid caught in the narrow gap between the two 
spheres is able to escape. The rest accumulates in a toroidal pocket, or bubble, of radius 
ri, oc rfl 2 that forms at the meniscus. This bubble is connected to a thin neck of width 
r n (x r 2 m . We have performed extensive simulations for the simplest case of equal viscosity 
fluids which confirm these scaling laws. 

In section IV we examine the inner "bubble" solution in greater detail. The bubble is 
joined to the neck by a short region of very large curvature on a lengthscale that appears 
numerically to be proportional to (r&rVj) 1 ' 2 and thus tends to a corner as r n /r{, — > (i.e. 
r m — > 0). Though the curvature of this corner is much greater than that of the bubble, both 
contribute at the same order to the leading-order motion of the meniscus, which can be 
thought of as simply due to a ring force of strength 2j smeared over a lengthscale rv In the 
final section we discuss the case of arbitrary viscosity ratios and mention related problems, 
namely, the effect of arbitrary initial shapes and the scaling at zero outer viscosity. 

II. FROM THREE TO TWO DIMENSIONS 

For simplicity we consider two initially spherical drops of equal radii R, as shown in 
figure |l|. Simple extensions to the cases of unequal radii and non-spherical shapes will be 
described in Section V. We denote the viscosity of the drops by r\ and that of the outer fluid 
by A -1 ^. As we have noted, the dynamics immediately after coalescence is described by the 
Stokes equations. Since these equations are linear, the velocity field can be expressed as an 
integral of the driving surface forces jkii, where n is the normal directed into the outer fluid 
and k = V.n is the curvature of the interface. We make the velocity dimensionless with 
respect to 7/77, all lengths with respect to R, and times with respect to the corresponding 
timescale r = -R77/7. 
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Calculation of the evolution of the interface S(t) requires only the interfacial velocity, 
which is given by the integral equation (Rallison & Acrivos 1978) 

(1 + A- 1 



u(xi) = - / Kj.ndcr 2 + (1 - A -1 ) / u.K.ndo-2, (1) 

Js(t) ' Js(t) 



where 



J < r) = SF 



I rr 

- + ^3 



3 rrr 

K(r) = -- — , r = Xl -x 2 , (2) 



rf<72 denotes a surface area element at position x 2 , and Xi,X2 both lie on S(t). The first 
term on the right-hand side of ([I]) represents the driving by the surface forces, while the 
second accounts for the difference in viscosity between the fluids. The problem is closed by 
requiring that any material marker £ with position x x on the surface moves according to 

9 iXl (0 = u( Xl ). (3) 

Equation ([]]) and the identity 

'j.nda = 0, (4) 



s 

(which is a consequence of the incompressibility condition V.J = 0) show that there would 
be no flow if k were constant over S. It follows that, in the early stages of coalescence when 
r m < 1, the flow is driven by the small region around the meniscus where k is not close to 
its initial constant value of 2. This key observation motivates an analysis based on splitting 
the interface into two regions. 

We use cylindrical polar coordinates (r, z) with origin at the junction between the drops. 
Away from the region of coalescence, the surface is essentially undisturbed and thus has the 
form h(z) = {2z) 1 ^ 2 and h(z) = (—2z) l l 2 for h 1 in z > and z < 0, respectively. The 
width of the gap between the spheres is given by 

w = r 2 (r m <r<l). (5) 

and, since dw/dr 1, the interfaces on either side of the gap are nearly parallel. 



The solution for this outer region has to be matched with an inner solution on the scale 
r = r m of the bridge or meniscus where the two drops are joined. The inner solution has 
a region of very high curvature, which provides the dominant contribution to the velocity. 
To a first approximation, this region can be represented as a ring of radius r m and small 
width A connected to two asymptotically straight interfaces each pulling outward with unit 
tension. The resultant effect is that of a radially directed ring force with strength 2 per unit 
length of the ring applied over a width A. 

To find the velocity field generated by this ring, we try integrating over a circular line C 
of forces f(r) = 2e r . Considering, for the moment, only the simple case A = 1, for which u 
can be computed directly from ([!]), we have 



u(xi) = — 

V 1J 871JC 



f(x 2 ) (Xj - x 2 )(f(x 2 ).(xi - x 2 )) 
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da 2 . (6) 



|Xi-X 2 | |Xi-X 2 

From this representation it is evident that the force distribution cannot be represented by a 
line everywhere as the first term in the integral (|6]) would lead to a logarithmically infinite 
value of the velocity. In the neighborhood of r m , when |x 2 — xi| = 0(A), one must account 
for the fact that the force is distributed over a length scale of size A. The logarithmically 
dominant part of the integral, which comes from A <C |x 2 — xi| <C r m , gives a radially 
directed flow 

u(r m ) = --!- In (—) e r . (7) 

2tt \r m J 

Since the curvature of the ring is not apparent at leading order in the region A <C 
|x 2 — xi| C r m that dominates ([?]), we may equivalently consider the corresponding two- 
dimensional problem, in which coordinates (x,y) take the place of (z,r). In that case (two 
parallel cylindrical drops connected along a narrow band of width 2r m ) the high- curvature 
meniscus is represented by two straight lines a distance 2r m apart. Since the forces 2j on the 
lines pull in opposite directions, they cancel on scales much greater than their distance apart, 
the integral (|B|) is cut off on the scale r m , and (0) again results. Because of this asymptotic 
equivalence of the two-dimensional and axisymmetric problems, we will mostly consider the 



two-dimensional problem from now on, which is simpler numerically. The two-dimensional 
forms of the kernels J and K can be derived by integrating (|2|) along the third dimension to 
obtain 



-Imr H — - 
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K(r) = r = Xl - x 2 , (8) 



7r r 4 



and the surface integral (|l|) is now along the perimeter of the two-dimensional drops. 

Not only does (0) give the leading-order velocity in both two and three dimensions, but 
it also holds for all viscosity ratios A. A summary of the argument is as follows. It is clear 
that the early flow can always be thought of as driven by ring or line forces of strength 2e r or 
2e r Since Stokes flow has no inertia, this force is transmitted unaltered across any surface 
enclosing the bubble and, since the width of the gap between the spheres is asymptotically 
negligible, the force must be supported by the internal fluid, and the external fluid in the 
gap makes little difference. Now the logarithmically large velocity of a slender body moving 
under a given force does not depend on the viscosity of the body, as can be seen explicitly 
in the solution for the motion of a cylinder of one fluid through another fluid (Lister & Kerr 
1989). Applying this result to the bubble at the meniscus, the finite viscosity of the external 
fluid in the bubble also makes little difference and (|7|) is correct at leading order for A ^ 1 
(though the higher-order corrections do depend on A). 

To evaluate the velocity from ([?]), however, it is necessary to determine the scale A over 
which the force is distributed at the meniscus and here the viscosity ratio does play a role. 
In the following section we will show that A ocr^, where a = 3/2 for finite A and a = 3 for 
the special case A = oo (no outer fluid). By integrating (0), we find that 

(en — 1) , 

r m {t)~- n J tlnt. (9) 

Recalling that time is measured in units of rjR/j, we see that the estimate based on di- 
mensional analysis alone, r m cx £7/77 is not quite correct, and in fact requires a logarithmic 
correction. 
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III. ASYMPTOTIC SHAPE OF THE MENISCUS 



Using the equivalence of the two- and three-dimensional problems, we now study the 
coalescence of two viscous circular cylinders in more detail. Figure |^ compares Hopper's 
exact solution for A = oo with a numerical simulation for A = 1 in which the initial condition 
is that of a cusp with r m = 10~ 3 , smoothed on the scale of r^. (The numerical method is 
described in more detail below.) The shape of the meniscus for A = oo is the tip of a near 
cusp (Fig. |2|a), while the shape for A = 1 is observed to be quite different for most of the 
evolution (Fig. ^|b): the external fluid is collected in a small bubble at the meniscus, making 
the lengthscale A of the local solution much larger than in the absence of an external fluid. 
Only in the last stages of merging does the fluid caught inside the bubble escape and the 
results look qualitatively more like Hopper's solution (Fig. 0c). 



A. Analysis for an inviscid exterior 

The existence of an exact two-dimensional solution for the special case A = oo of an 
inviscid or absent external fluid (Hopper 1993a,b; Richardson 1992) allows us to test the 
general ideas of section II. Asymptotic expansion of this solution near the meniscus shows 
that 



h(x) 



2^ + V(l r m) +( 2X 



1/2 



(x<l), (10) 



from which we deduce that the highly curved region is of size A ~ r^, and the curvature 
k scales like A -1 ~ r m 3 . This result is somewhat surprising since A is much smaller than 
the gap width w ~ estimated from the spherical shape of the outer solution and, at 
present, we do not have an asymptotic argument for the appearance of the small scale r 3 ^. 
(By contrast, with a viscous outer fluid A ~ r 3 / 2 is larger than the gap width w, and this 
scale can be understood from mass conservation as discussed in section III B.) 
Inserting A ~ r^ into (^), we obtain 
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v ~ lnr m (11) 

7T 

for the velocity (v = r m ) at the meniscus, which agrees with the asymptotic result given by 
Hopper (1993a); our earlier asymptotic analysis now allows the extension of this result to 
three dimensions. 



B. Analysis for a viscous exterior 

For the case A = 1 we begin by discussing the structure of the local solution close to 
the meniscus, which is shown in figure ||]. It consists of a "bubble" of outer fluid of radius 
A = r&, which is connected to a thin neck of width r n . The neck matches onto the static 
outer solution, so r n must scale like the gap width w ~ r^. The area of the original gap up 
to r m is O(rj^) so that, if the meniscus advances faster than fluid can escape from the gap, 
the bubble should contain a finite fraction of the gap fluid and hence cx rfl 2 . 

To examine this argument in more detail, we consider the velocity field generated by the 
large curvature of the meniscus. As we have already noted, the flow is driven by that part 
of the interface where the interfacial curvature is significantly different from 1 (or 2 in the 
spherical case), namely \r — r m \ = 0(r&). Using a multipole expansion of this forcing, we 
find the velocity field u(xi) at a distance |xi — x m | ^> r& from the center x m of the bubble 
to be 

u(x 1 ) = f.J(x 1 -x m ), (12) 

at leading order, where J(r) is given by @ and f is the total force exerted by the bubble. 
This force is the integral of — kii = d s t over the bubbles' surface (£&), where t is the tangent 
vector pointing in the direction of increasing arclength s: 

f = — / /mdcx = / d s tda = t 2 — ti. (13) 

In the present case — ti = t 2 = e y , so that the far- field velocity resulting from the forcing of 
the bubble and its image is 
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u(xi) = 2[J(xi - x m ) - J(xi + x m )].e y (14) 

We are interested in the flow in the neck, so we choose Xj = (0,yi) with y x — r m ^> r b . 
From @ and ( |T4"1) we find that the y-component of velocity in the neck is given by 



«n(yO = — 1b(?^— ( 15 ) 



2tt Wi + r 



The representation (FJ) breaks down when y\ — r m = 0(r&), since the higher-order terms in 
the multipole expansion become of comparable magnitude and u n (yi) crosses over to some 
function that depends on the detailed structure of the bubble. Since (|i~5D must match onto 
the velocity field in the bubble, we can write the velocity of the meniscus as 

1 , /2r 



v = v + ^ln(^\, (16) 
27r V r b J 

where the constant v$ comes from the detailed shape of the bubble, and for A = 1 is found 
numerically to be v = —0.077. (This equation is consistent with (|7|), but also includes a 
representation of the 0(1) contribution.) 

If rj <C i/i - r m <C r m then from ( jl5[) and (|T6|) 

»n(m) * ^ In {-^A « i- In (Z=) « v. (17) 



2tt \y x - r m J 2vr V r b 
Thus the fluid in the neck ahead of the bubble only starts to move at a speed comparable to 
that of the bubble when yi — r m = 0{r^) i.e. when the bubble has caught up with it. Thus 
all the fluid in the neck is collected into the advancing and growing bubble. Since the neck 
width r n scales like r^, the total area of neck collected scales like ~ r^, and thus 

n-r^ 2 . (18) 

Finally, combining (18 ) with ( [TBD or @, we find that 

v ~ -— lnr m , (19) 

47T 



which can be integrated to give 



r m ~-— Int. (20) 

47T 



This result differs by a factor of 4 from the A = oo limit. 
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C. Numerical tests of scaling 



To test the predictions of these scaling ideas, we have performed extensive two- 
dimensional simulations of drop coalescence using the boundary-integral method. The initial 
bridge has r m (0) = 10 -6 , which means the gap width is w = 10~ 12 initially. For simplicity, 
we only considered the viscosity matched case A = 1 so that the second term drops out from 
(H) and u can be computed directly from the surface forces. The interface was parameterized 
by arclength, and derivatives were evaluated using centered differences. The interface was 
advanced according to using an explicit second-order Runge-Kutta step. The difference 
between the result of time step At and two half-steps At/2 was used to control the time 
step. 

Improvement of the stability of our numerical method by making it implicit would be 
computationally very demanding for an integral operator. Instead, we use a scheme first 
proposed by Douglas & Dupont (1971). The equation of motion (^ for the position of the 
interface can be written at any given time as a linearization around the current interfacial 
position x : 

<9tx = A(x) ■ (x — xo) + constant (21) 

By writing this as 

d t 5x = (A - B)5x + B5x + constant (22) 

and treating the first part explicitly, but the second part implicitly, the scheme becomes 
unconditionally stable as long as |B| > |A|/2, where the matrix norm is defined to be 
the modulus of the largest eigenvalue (Douglas & Dupont 1971). In the present case, |A| 
scales like (Ax)~ l up to logarithmic corrections, where Ax is the minimum grid spacing. 
By choosing B to be a diffusion operator multiplied by the local grid spacing, one can make 
sure that the numerical method, although treating the integral operator explicitly, becomes 
unconditionally stable. Without the help of this trick the time steps required to integrate 
over the first decade and a half in r m would have been prohibitively small. 
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To achieve the necessary spatial resolution, local refinement of the mesh is crucial. The 
resolution near the meniscus was set by the inverse of the local curvature. Away from 
the inner solution, the grid spacing was allowed to taper off geometrically, with the spacing 
constrained to change by no more than 10% from one grid point to the next. As explained in 
more detail below, additional resolution was used in the transition region where the bubble 
merges into the neck. The maximum number of points used to represent one quadrant of 
the shape was about 900. Every few time steps a new grid was constructed using the current 
interface, and the interface was interpolated to the new grid. Thus there was no need to 
rearrange the grid points along the interface. 

In figure |] we show the scaling of the bubble radius r b and the minimum neck radius 
r n and observe that both follow the predicted power laws (0). A closer inspection shows 
that the slope of logr^ is slightly smaller than expected, which is because the scaling of the 
area of the bubble is almost the same as that of the neck. (The portion of the neck up 
to the point y = r m + ar b contributes an area A n « a^V^i ~ r^ 2 , which is only slightly 
smaller than Aj, ~ since the two exponents are close.) We confirmed numerically that 
the r m -dependence of the total area A = A b + A n has no significant deviation from r^, as 
expected. 

The scaling of the velocity at the meniscus (v) and in the neck (u n ) at position y\ = 
r m + art, is shown in figure |5] (with a = 20) and compared with (1/2tv) ln(2r m /r;,) , which 
is predicted by equations ( fL6|) and ( |PTD to have the same slope. The theory gives v — u n = 
constant = \na/(2n) + t> which is found numerically to be very close to 0.4 and hence 
v = —0.077. All three curves in figure || should have a slope In 10/47T when plotted against 
log 10 (r m ). The noticeable deviation comes from the fact that one is effectively taking the 
difference between \og 10 rfl 2 and log 10 r m , so non-asymptotic effects in r b , still present on 
these small scales, become more pronounced. 

It has been implicit in the previous arguments that the inner solution, consisting of a 
bubble connected to a thin neck, has reached its asymptotic form: in the frame of reference of 
the advancing bubble tip, and rescaled by the bubble radius, the shape should be stationary. 
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Figure || showing both the local interface profile and the curvature at two values of r m one 
decade apart, reveals that there is, in fact, a slow variation in part of the local profile. This 
slow variation is seen as a positive second peak in the curvature at the point where the 
bubble meets the neck. To be able to resolve scales down to r m = 10~ 6 , additional grid 
points were inserted at the position of the second peak, where the grid spacing was based 
on the width of the peak. As r m — > 0, which corresponds to going back in time, the second 
curvature peak increases and also gets narrower, as its integral must be finite to yield a finite 
change of slope. 

To obtain more information about the asymptotic shape of the inner solution and the 
growth of the secondary peak in curvature, it is useful to consider the inner solution as a 
separate problem. This analysis is given in the next section. 

IV. BUBBLE ON A NECK 

Here we study the local solution close to the meniscus, which consists of a bubble of 
radius r& connected to a thin neck of width r n (see inset to figure [7]). Asymptotically, 
the curvature of the neck is very small compared with , so the neck can effectively be 
considered as an infinitely long channel of uniform width. From now on, all lengths will be 
measured in units of r?,, so the radius of the bubble is normalized to unity, and the radius 
of the neck asymptotes to some small number e ~ r n /rb. The solution we are interested in, 
which corresponds to the asymptotic structure of the main solution described in Section III, 
is such that the interfacial shape g(y, t) is advected at a constant speed v c without changing 
its shape 

g(y,t) = G(y-v c t). (23) 

The physical meaning of this statement is that in the original problem the local solution 
relaxes to a quasi-steady shape on a much shorter timescale than the position r m of the 
meniscus is changing. The velocity field of the local steadily translating shape ( |2"3"| ) must 
satisfy 
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(uy - v c )G' = u x . (24) 

The boundary condition is that g(y) approaches e as y — ► oo. The components of the velocity 
field u x and u y follow from the integral equation ([1]) as usual. Instead of solving the system 
([!]) and (^4]) directly, we found it most convenient to evolve two bubbles attached to a very 
long, straight neck of radius e^n until a stationary shape is established, as shown in the 
inset to figure |7|. The tension in the neck is responsible for pulling the bubble along. The 
neck shortens during the relaxation and the radius increases to a value e, which then only 
changes very slowly by the time a stationary shape is reached. 

Since the radius of the bubble and the radius of the neck r n are very different, one 
might think that the curvature distribution within the bubble is independent of the neck 
radius e. However, the limit e — > turns out to be singular as an increasingly pronounced 
peak of positive curvature appears at the junction between bubble and neck, as demonstrated 
in figure |6|. 

First, in figure [7] we compare the curvature distribution as given in figure |6] for r m = 10~ 3 ' 5 
with that of the stationary problem, equations (|TJ) and (pi"D, with e = r n /r b . We choose e init 
such that the neck has the appropriate width e by the time a stationary shape is reached. 
The excellent agreement shows that the flow close to the meniscus is completely equivalent to 
the translating bubble, which is of course a much simpler problem. Hence the inner solution 
of the coalescence problem can be understood completely in terms of the translating bubble. 

Figure |8|(a) shows a sequence of bubble shapes for increasingly small values of e. The 
overall shape of the bubble does not depend very much on e, and it looks as if the shapes are 
almost the same for the two smallest values. However there is an increasingly sharp "corner" 
at the point where the neck meets the bubble. This result is most evident from a plot of 
the curvature for the same values of e; see figure |8|(b). While the curvature of the bubble 
is negative, the corner at the junction between the bubble and the neck corresponds to a 
growing peak of positive curvature. In figure |S](b) we also include a plot of the maximum of 
this peak as a function of e. The data is suggestive of k' 1 ~ e 1 / 2 though it is hard to make 
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a definite statement. If this is so then the lengthscale of the corner is the geometric mean 
of the scales of the bubble and the neck, which would be r^ 4 in the coalescence problem. 
From ( |i3|) it is evident that the total force exerted by the peak is equal to the change in 
slope. As can be seen from figure || (a) the change is constant to a good approximation so 
that, of the total force 2e y exerted by the bubble, roughly 15% is exerted by the corner and 
85% by the rest of the bubble. 

V. DISCUSSION AND OUTLOOK 

As we have demonstrated, even the simplest case of viscosity-matched fluids represents a 
problem of enormous complexity, in which there are features on at least the lengthscales r 2 ^, 
r rl 2 i r m an d 1. Therefore, we confined ourselves to computing the leading-order asymptotics, 
which we note are only logarithmically dominant. In all likelihood, quantities like the 
bubble radius J*&(t) contain additional logarithmic terms, whose calculation require a better 
knowledge of the matched asymptotics of the problem. Further complications arise because 
the inner solution is itself singular, with a corner of lengthscale which we estimate to be 
0(r^ 4 ). Important goals for future work would be confirmation of the inner scalings, formal 
asymptotic matching of the different scales, and to go beyond the leading-order problem. 

Another important problem is the generalization of our calculations to arbitrary viscosity 
ratios A. A major obstacle to developing a quantitative theory for general A is that numerical 
simulations are much more difficult. When A ^ 1 (jlD is a second-kind integral equation for 
u, which requires an order N 3 effort to solve instead of N 2 . Moreover, for small r m the 
matrix associated with the second integral becomes singular since the local solution for the 
velocity field is close to a uniform translation for which the kernel has a zero eigenvalue. We 
have not yet found a way to treat this singularity sufficiently well to go beyond r m -values of 
lO" 2 for A ^ 1. 

We have argued already that we expect the leading-order behaviour of r m to be given 
by (|7|) for any A since the net force from the meniscus is supported asymptotically by the 
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internal fluid. We also believe that the scaling A ~ r 3 / 2 will hold as r m — > for any finite A 
since, even for large A (small external viscosity), the pressure drop along a narrow channel 
is large and so the outer fluid is not able to escape from the gap and is caught in a bubble. 
Expressed differently, the positive curvature at the bubble corner is unable to pull the walls 
of the channel apart since this would require significant motion along the channel. It is 
intriguing that Hopper's exact solution shows that the situation is very different for A = oo 
when there is no external viscous resistance to overcome and no bubble forms (figure 2a). 
Thus the limits A — > oo and r rn — > do not commute, reflecting the fact that one must be 
careful in assuming a zero-stress condition in situations involving narrow cusps, as even a 
very small external viscosity can be significant. 

We can provide a physical estimate for the scale below which a bubble forms for A ^> 1. 
Fluid motion in the narrow gap may be treated with the lubrication approximation, whence 
r]u n /(Xr 2 ) ~ Ap/r m and Ap « 7/r b gives the estimate u n ~ (\^ / r]){r m / R^l 2 . Therefore, a 
bubble can only be expected to form when the meniscus motion v ~ j/r] > u n , which occurs 
on an approximate lengthscale r m j R < A~ 2 / 3 . 

We have been considering the simplest case of equal spheres brought into contact. The 
case of unequal spheres of radii R and R/5, with 5 < 1 is a straightforward generalization 
in which the gap thickness (|5|) is simply replaced by 

w = r 2 (l + 5)/2 (r m <r«l). (25) 

Similar considerations show that the initial evolution of axisymmetric drops brought into 
contact along their symmetry axes depends only on the local curvature at the initial contact. 
An more interesting variation is that of general initial shapes for which the locus of high 
curvature (in three dimensions) near contact no longer forms a circle, but is a more general 
closed curve. To leading order, this curve is convected with a logarithmically large velocity 
field, pointing in the direction normal to the curve. 

Recent research (Nikolayev et al. 1996; Bonnecaze et al. 1998) has suggested that the 
rate of coalescence in emulsions can be greatly enhanced by the flow generated by individual 
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coalescence events. Our analysis suggests that an appropriate model for the far-field velocity 
of a single coalescence event is the Stokes flow driven by an expanding ring force of radius 
r m and strength 27 per unit length. This dipolar flow, which may be obtained by solving 

r/V 2 u - Vp + 27<5(r - r m )5(z)e r = (26) 



using Hankel transforms [22], is 



u r {r,z) = ^ r J 1 (kr)J 1 (kr m )(l-kz)e- kz dk 



2r] Jo 

u z {r,z) = r J {kr)J 1 {kr m )kze- kz dk. (27a) 

2r\ Jo 

p(r,z) = 2i]u z (r,z)/z, 

where the J n are Bessel functions. These equations can be approximated in the limit r m <C 
(r 2 + z 2 ) 1 ^ 2 to obtain the axisymmetric dipole 

which is shown in figure |^ superimposed on the outlines of two spherical bubbles. The 
coalescence-induced radially directed force drives a flow towards the two drops over an angle 
w 127°. The flow Q2"8|) might be used for simplicity in models of multiple coalescence. 

All motions described so far begin with a local point contact and it is worth considering 
how this contact might be achieved. The near-contact squeezing motion generated when two 
drops (or a drop and a plane) are in relative motion can be analyzed using the lubrication 
approximation (e.g. Jones & Wilson 1978; Yiantsos & Davis 1991). Owing to the large 
pressures accompanying flow along the narrow gap, the surface tends to deform in the 
narrow gap. In particular, when two equal size drops are squeezed together with a force F 
on each drop, then a dimple tends to form when h (t) < R/X 2 , where ho(t) is the minimum 
gap spacing and R the radius of curvature of the undeformed drop. The magnitude of the 
deformation becomes the same order of magnitude as the gap height and the dimple has 

1 /2 

a radial scale (FR/ r y) 1 . Away from the gap the drop is nearly spherical so long as, for 
A = 0(1), the effective capillary number is small, 0(Fh^ f^fR?/ 2 ) 1. An implication of 



these analyses is that contact is very likely to occur along a rim, or at least at an off-axis 

1 /2 

position, with an initial radius of the bridge of order (FR/j) . The dynamics would then 
follow the results shown in Section III. 

From an experimental point of view, it is probably not relevant to investigate smaller 
(dimensionless) r m than 10~ 5 since the gap width just ahead of the bubble is proportional to 
r 2 ml which is then of microscopic size for reasonable values of R. Very small inhomogeneities 
in the fluid or van der Waals attractions will cause the two interfaces to reconnect, and to 
create an instability that breaks the azimuthal symmetry we have assumed in the three- 
dimensional problem. Moreover, the bubble actually forms a structure that resembles a long 
thin torus in three dimensions, and is thus prone to a Rayleigh capillary instability, which 
grows on a short timescale proportional to and is potentially dangerous. On the other 
hand, there are stabilizing effects since the bubble is also convected (Brenner, Shi & Nagel 
1994), and a careful nonlinear stability analysis has to be done to determine the scale where 
the stability is first expected to occur. This question is experimentally relevant because it 
sets the size of small bubbles of the outer fluid that may be observed after coalescence. 

We have pointed out that coalescence is initially described by the Stokes equations. If 
the viscosity of the fluid is small, this is true only for the early stages of coalescence, until 
the Reynolds number is of order one, which happens when 

r m ^i, £ = v 2 /(p 7 ). (29) 

For water, I = 1.4 x 10~ 6 cm , so it is a very relevant question to go beyond the Stokes approx- 
imation. After passing the transition region (29) we expect the dynamics to be described by 
the Euler equations. Assuming that the scale of the local solution at the meniscus is set by 
the gap width alone, the interfacial stress driving this motion is approximately ^/(r^/R), 
which is to be the same magnitude as pv 2 . Hence, with v = r m , we find 

r m oc (^V 2 , (30) 

which corresponds to v oc t^ 1 ^ 2 . The geometrical part of both problems is similar to the 
Stokes case, but the coupling between pressure and velocity makes the relationship between 
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the surface shape and the velocity different, so an alternative treatment is needed to predict 
the numerical coefficient in front of the power law for r m . 

Clearly, the possibility of finite Reynolds number, arbitrary surface shape, finite velocity 
of approach, and the inclusion of another fluid outside the drops lends a tremendous richness 
to the class of singularities studied here. 
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FIGURES 

FIG. 1. The surface profile h(z) produced by two coalescing drops of radius R. The origin of 
the axis of symmetry z = lies at the initial point of contact. The bridge joining the two spheres 
has radius r m and width A. 

FIG. 2. A closeup of the point of contact during coalescence of two identical cylinders for 
the two cases of no outer fluid (A = oo) and two fluids of equal viscosity (A = 1). (a) is Hopper's 
solution for r m = 10~ 3 , 10~ 2 ' 5 , 10~ 2 , and 10 -1-5 . (b) a numerical simulation of the viscosity matched 
case that shows fluid collecting in a bubble at the meniscus. Note that the two axes are scaled 
differently, so the bubble is almost circular. For large values of r m , as shown in (c), the fluid finally 
escapes from the bubble, and the width of the meniscus is closer to the value of the gap width 
0(rl). 

FIG. 3. The structure of the local solution close to the meniscus. It resembles a bubble 
connected to a thin neck. The radius of the bubble is r& and the minimum radius of the neck is r n . 
The distance from the origin to the front of the bubble is r m , which is not drawn to scale here. 

FIG. 4. Scaling of the bubble radius r& and the neck radius r n as function of r m . 

FIG. 5. Scaling of the velocity v at the tip and the velocity u n at a position y = r m + 20r& 
in the neck. There is a constant difference of 0.4 between the two. Both agree very well with 
the scaling of (l/27r) ln(2r m /r&) as predicted by theory. For comparison, we also give a slope of 
(lnl0)/47r = 0.183. 

FIG. 6. The local solution and its curvature distribution in coordinates rescaled by r&. As r m 
gets smaller, a sharp peak develops at the junction between the bubble and the neck. 
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FIG. 7. A comparison of the curvature distribution of the local solution for r m = 1CT 35 and 
the stationary "bubble on a neck". The peak height of the positive curvature has been used to 
match the two. In an inset, we show the initial condition used to compute the stationary shape of 
the translating bubble. 

FIG. 8. (a) The stationary state of a translating bubble, pulled by a thin neck for different 
values of — log 10 e = 1.8,2.8,3.7, and 4.9. (b) Blowup of the curvature distribution. Although 
the surface shape seems to have converged, the maximum negative curvature in the junction still 
increases. The insets show the scaling of the maximum curvature n max and that of n max * 1^1/2, 
where Wi/ 2 is the half-width measured in arclength. 

FIG. 9. The velocity field generated by an expanding ring of forces in the limit of the radius 
of the ring r m going to zero. Note that the velocity field is pointing outward in the direction of 
the expanding ring, but is inwardly directed over much of the flow domain. 
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Figure 2c: 
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Figure 5: 
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Figure 8a: 
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Figure 8b: 
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